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Abstract 

We show that in theories where the lowest energy excitations are not 
quasiparticles but they form a continuum, the shear viscosity to entropy 
density ratio goes to zero as the temperature goes to zero. In these the- 
ories therefore there is no lower bound for the shear viscosity to entropy 
density ratio, in contrast to the predictions coming from the AdS/CFT 
correspondence. 
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1 Introduction 

Experimental evidences from RHIC suggest that the plasma formed after heavy 
ion collisions is very close to a perfect fluid. The measured large value of v% 
[2 [5] indicates a streaming with very small viscosity [31 H] : fits to the values 
suggest an upper bound r\j s < 0.16 [5]. 

There is a tremendous effort in the recent literature to give an account for 
such small values of viscosity from the theory side. Perturbation theory a priori 
performs an expansion around the ideal gas limit, where 77 = 00, and so the weak 
coupling expansion usually yields large rj/s ratios proportional to l/g 4 lnl/g 
[HI HI M ■ Still there is hope to access experimentally acceptable ratios using 
higher order processes in Boltzmann equations [101 111) . There exist analytic 
methods beyond weak coupling perturbation theory, like large Nt expansion 
|12j or strong coupling expansion in pure Yang-Mills theories |13j . 

One can use Monte Carlo simulations to calculate energy-momentum tensor 
correlators, which is related to the viscosity via an integral equation. From this 
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one can extract the viscosity for pure Yang-Mills theories [131 [23 HH] • These 
results are, however, not fully reliable yet [T7] . 

A method proved to be very useful in studying the viscosity to entropy ratio 
is based on the AdS/CFT correspondence: in J\f = 4 super Yang-Mills theory 
at strong t'Hooft coupling A and at large N c one finds r\j s = 1/4tt [18 . In lot 
of other gravity models we find the same result [TS] . One can also calculate the 
1/A corrections 20 , which turned out to be positive. It is very natural then 
to expect, that there is a lower bound to the viscosity to entropy density ratio 
in the conformal case. One can argue that this may be true for other models, 
too [21j : since the viscosity over entropy ratio parametrically is proportional 
to Et, where E is the quasiparticle energy, r is its lifetime, we expect that 
Heisenberg uncertainty relation constrains this ratio. The value of the lower 
bound is expected to be given by the conformal case, l/4n. 

There are, however, caveats in the argumentation, in two aspects. One is the 
gravity side: in higher derivative gravity models the ratio can be smaller than 
1/47T. In the construction of Ref. [52] the lower bound in the unitary regime 
is 16/25 x 1/4-7T [23]. There are other string theory constructions, where one 
can violate the 1 /An lower bound [24] . Recently there are several gravitational 
models, where the value of the lower bound is challenged further J55[[2n]. This 
may mean, that in a generic higher derivative gravity model there is no lower 
bound at all. 

The other caveat comes from the particle physics side: large number of 
particles can lead to a large Gibbs mixing entropy, while the viscosity is constant 
[27l EU [29[ [12] , and so the rj/s ratio approaches zero. In [30. it is shown that the 
large number of species can be in the same quantum channel as excited states. 

In quantum field theories the density of states generally contains a continuous 
part, which can be considered as the generalization of the large number of 
excited states. In these models under certain conditions we may expect the 
violation of the rj/s bound. This was observed in strongly coupled pure Yang- 
Mills theories [T3J: when going to zero temperature the viscosity to entropy 
density ratio approached zero. From the details of the calculation one finds 
that this behavior is the consequence of that that in the spectrum there are no 
quasiparticles, only a continuum above a mass gap. 

In this paper we show that this behavior is a common feature of a large class 
of relativistic quantum field theories, where the lowest lying excitations of the 
energy spectrum form a continuum instead of a discrete set (quasiparticles). As 
opposed to systems with quasiparticle excitations (Bose gases or liquids), these 
bosonic systems are non-Bose gases (liquids). In these models the rj/s ratio 
approaches zero when the temperature goes to zero. In fact the argumentation 
is based on kinematics, and therefore it is even more generic: we should find 
this behavior for all transport coefficients. This result may suggest that these 
theories have, if they have at all, a nonlocal dual gravity model. 

Before the actual calculation we note that formally all relativistic theories 
at zero chemical potential are superfluids, since the spectral function must be 
zero below the light cone. It does not mean that the viscosity goes to zero at 
zero temperature: in fact also in superfluid He the shear viscosity is divergent 
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[3Tj . As we will show the T — ► limit of rj/s can be very different for different 
models: it can be infinite (massive quasiparticles), finite (conformal case) or 
zero (non-Bose liquids). 

2 General setup 

To start we consider a self-adjoint operator A(x), and study its correlator at 
low temperatures. We use the standard technique [35] to write it with the help 
of matrix elements of A(Q). To this end we choose a basis in the Hilbert space 
and write at finite temperature 



C A (x) = ([A(x),A(0)}} = 

Tr(n \e- pH A(x)\ m )\n) 



Z^ 

n,m - 

Tr ( to \e' l3H A(0) \n)(n \A(x) | m } 



(1) 



where Z = Y^neN Tr e~P H . For the states we assume that they are represented 
by the four- momentum of the state p = (po>p)j an d other quantum numbers 
Q, so we have \n) — \p, Q). Note, that p and p are not connected in general 
by a dispersion relation, not yet the possible values of po form a discrete set. 
To describe the generic situation we introduce the density of states of a given 
quantum channel Qq(p), with the definition 

X>|..>) = W^ eQ (p)< Pl G|...|p J e). (2) 

We will also assume that the normalization of the states are done for unit 
volume - then all calculated quantities are densities. We emphasize here that 
the volume of the system is infinite, only the normalization volume of the states 
is fixed! 

We can write A{x) = e ipl * x >* A{Q)e~ lP ^ x ^ with the use of the generator of 
space-time translation, and so (q,Q\A(x)\p,V) — e~ l ( p ~ q ) x ( q, Q \A(0)\ p, V ). 
We also perform Fourier transformation, and find finally 



x (1 - e-»") (2*fHk + q - p) !(„, Q |A(0)| p,P )| 2 . (3) 



This form demonstrates also that C(fc) > 0Vfc > 0, and C(—k) = —C(k). 

The transport coefficient corresponding to A, which is denoted here as 77^4, 
comes from the Kubo formula [55] : 

i] A = lim — . (4) 

fc ->o fc 
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It reads with the current representation 

E J^QQ(l)Mq)e^\(q,Q\A\ q ,V)\ 2 

(5) 



where we have omitted the zero argument of A. 

Now we go to low temperatures, (3 — > oo. The exp(— {3q ) factor forces to 
keep only those quantum channels, which contain the lowest energy levels while 
the matrix element is not zero there. For simplicity we assume, that this is true 
for a single quantum channel, but this is not crucial for the later discussion. 
Therefore in the followings we will suppress the notion of the quantum channel, 
and write 

U^r A 9 2 {<l)e-^\(q\A\q)\\ (6) 



VA ZJ (2tt)4 

Analogous result could be found for the viscosity in strong coupling expansion 
at small temperatures |13j . 

We should work out the low temperature expression for the entropy, too. 
The free energy is the logarithm of Z, which can be written, by explicitely 
separating the vacuum contribution: 

z-i= E <-l^l-> = E f£fre*to'- fi *- ^ 

At low temperatures the correction to 1 is small, and so we can approximate 
ln(l + 5Z) = 8Z to leading order. We again assume that the most important 
contributions come from a single quantum channel, and suppress the index Q 
in the followings. It could also happen, that the relevant quantum channel for 
entropy is different from the one for r]A, but we will consider only those transport 
coefficients, where this is not the case. Transport coefficients coming from the 
commutator of the energy- momentum tensor (like the shear viscosity), are in 
this class, since all excitations necessarily contribute to the energy-momentum 
tensor. Then at low temperatures the free energy density and the entropy 
density read: 

At low temperatures, therefore, the ratio of tja and the entropy density reads 



dT J (2tt) 

where we have taken into account that at low temperatures Zal. 

In the following we make kinematical considerations to estimate the temper- 
ature dependence of the r\j s ratio. There we should assume something about 
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the matrix element |(g|^4|g)| . Since we calculate transport coefficients, we 
may assume that A is a current, and therefore its matrix element is zero, when 
there is no gradient of the corresponding density. Therefore we assume that 

|( g L%)| 2 = q 2 A (10) 

The coefficient A can be temperature dependent, but at small temperatures we 
may approximate it with its zero temperature limit. One should note that for 
the shear viscosity the current is a tensor, because we describe the current of 
the energy-momentum four vector. But it does not alter the statement that the 
matrix element remains to be proportional to q 2 where q describes the change 
of the momentum in the perpendicular direction. 
So we can write: 

- = ^7^ • (") 

s d m f 6rq , . Rn 

^ T /(4^ )e 
We will consider two different types of systems: one with a quasiparticlc behav- 
ior at the lowest energies, the other with a threshold behavior, and examine the 
tia/s ratio. 



3 Quasiparticle case 

Here we assume a Breit-Wigner-type distribution of the lowest energy eigenval- 
ues: 

2T 

M = (»-„,)■ + 1- <12 » 

where e q is the dispersion relation of the quasiparticlc, T is the quasiparticle 
width. Here we will assume a relativistic dispersion relation e 2 = q 2 + m 2 . If 
r — > this form approaches 2Tr5(qo — e q ). We assume that we are in the "small 
width" regime, which means that always, where we can, we should send r — > 0. 

This policy suggests that we should treat the product of two distribution 
function in the sense of the Fermi's Golden Rule, ie. we treat one of them as if 
it was a Dirac-delta: 

2 

q (q) ~ q(qo = £q)Q(q) = p QW- (13) 

Thereafter the integrals both in the numerator and in the denumerator contain 
one spectral function, which can now be approximated by the delta function. 
The free energy density reads 

oo 

/— r/^t «(«)«-*- -st/*^.-*-- 



oo 

= -— J dzzVz 2 -le- fjmz - -^([3m) 2 K 2 (pm), (14) 
l 
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where K2 is the modified Bessel function of the second kind. At zero mass 
(conformal case) the relevant limit is (5m — > 0, then we find 

_ T 4 4T 3 

J|conf — T, ^|conf — T~ • \*-^} 

TT 7T 

In the massive case at temperatures T <C m we find in the leading order: 



For the transport coefficient we find with help of (fT3|) 

■» - • 4a / 1? ^W'-" ^/l? = 

00 

• 4 " i ' [ dzz( Z 2 -lfV"^K 3 (W, (17) 



TrTr 2 7 v ; Ttt 2 
1 

where -K3 is the modified Bessel function of the third kind. In the conformal 
case: 

24AT 4 

?M| con f= r 2 ■ (18) 
In the massive case, at temperatures T <C m we find in the leading order: 

?M|m " YT (2tt) 3 / 2 ■ 1 ' 

Therefore the tja/s ratio reads in both cases 

5* = <ul (20) 
s I 

As we can see, the desired ratio is proportional to T/Y. This is the man- 
ifestation of the qualitative argumentation of Ref. [21] : since T is the kinetic 
energy and \/Y is the lifetime, therefore t]a/s ~ E^t. Therefore this ratio is 
constrained by the Heisenberg relation, and we expect to have a lower bound. 
In fact, in the conformal case Y must be proportional to the temperature (there 
is no other scale), and the ratio is constant. In the massive case, on the other 
hand, the width is exponentially small. This is because to form a width, the 
particle must scatter on thermal states, but their abundance is ~ e~ M ^ T , where 
M is the energy of the lowest lying scattering state. Therefore 



s 



VA 

~ const., — 

conformal 



Te M/T i^O ^ (21) 



G 



4 Non-quasiparticle case 



The other big class of theories is when the low energy spectrum cannot be de- 
scribed by a narrow quasiparticlc peak, but instead there is a continuous density 
of state. In this case it is worth to stay in the four-dimensional formalism, and 
write the integral in 4D polar-coordinates. We should take into account that, 
since the energy spectrum has a lower bound (for stability) which can be chosen 
to be zero, the qo integral is restricted to positive values. Moreover the integrand 
is spatially symmetric. Then we can use: 



/ -> i/ dqq3 J ^ sinhV (22) 

o o 

The expression for the free energy can be rewritten as 

oo oo oo 

/ = -^3 J dqq 3 g(q) J dq sinh 2 r,e~^^ = ~ J dqq 2 g(q)Ki(/3q). (23) 



For the transport coefficient we obtain 

oo oo oo 

VA = ^1 dqq 5 g 2 (q) J d V sinh 4 rye^"^ " = J dqq 5 g 2 ; (q)K 2 (/3q). 



(24) 

At T — > limit we will use the zero temperature limit of the spectral density 
which is Lorentz-invariant. Since the modified Besscl functions all decrease 
as K(x) — ► (7r/2x) 1 / 2 e _:E , the (3q argument in the above expressions enhances 
the lowest momentum part of the spectral function, ie. the threshold region. 
Sufficiently close to the threshold the spectral density can be approximated by 
a power law: so we will assume that the spectral density has the form in the 
relevant regime as 

g(q) = CG(q-M)(q-M) w , (25) 

where M is the threshold value and C is some constant. 

If there is no mass gap, then the expressions for / and r\A can be exactly 
performed, yielding for the entropy density 



_ C(5 + 

4tt 3 " \~2 



(26) 



and 

n A = ^2^T 5+2w 4 1+w T(l + w)T{3 + w). (27) 
The transport coefficient to entropy ratio is therefore 

^ = a w ACT 1+w (28) 
s 
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where a w = 3 2 W T(l + w)T(3 + w)/[(5 + w)T((5 + w)/2)T((3 + w)/2)]. 

If M 7^ then at T <C M we can use the asymptotic form for the modified 
Bessel functions, and obtain 

8 = ^1 CT ^J V T -+3/2 M 5/2 e -/ 3 M ! (2g) 



and 

fir 3^C* 2 r(2w + 1) 

VA - 



T 5/2+2w M 5/2 e -f3M (gg^ 



The transport coefficient over entropy ratio is therefore 

24. = a w ACT 1+w , (31) 
s 

where a w = ST(1 + 2w)/T(l + w). 

So finally from both the gapped and gapless case we obtained the same 
result: 

rjA _ Tl+W t-* (32) 
s 

for any realistic threshold behavior (w > —1). 

We remark here that the result could be guessed by dimensional argumenta- 
tion. The coefficient C must be dimensionfull: g is of dimension 1/E, q, M ~ E 
therfore C ~ E^ 1 ^ w . Since ija contains C 2 , the entropy density C, their ratio 
is ~ C. But the ratio is dimensionless, so a quantity with dimension E 1+w must 
appear. If M = the only candidate is T, so we must have T 1+w . If M ^ 0, in 
principle there could be also some M factor, but it turns out from the concrete 
calculation, that the M factors drop out. 



5 Conclusion 

For conclusion we recall that we studied the low temperature behavior of the 
ratio of some transport coefficient tja (where the notation is motivated by the 
shear viscosity) and the entropy density. The computation is based on exact 
formulae which are used near zero temperature. We studied two classes of 
models, both describing relativistic superfluids. 

The first is the class of the Bose-liquids, where the lowest energy excitation 
is a stable particle at zero temperature which becomes a narrow width quasi- 
particle at small temperatures. To this class belong a lot of particle physics 
models, for example the <f> 4 model, or the sigma model. There we have found 
that the tja/s ratio is proportional to T/T, in agreement with the qualitative 
expectations |21j . In the conformal limit this ratio is constant, in the massive 
case the ratio is diverging as T — » 0. In this class therefore there is a lower 
bound for the t\aI s: in the case of viscosity it is probably 1/4tt [T5] , 

The second class of the studied models is the non-Bose-liquids: here the 
lowest lying excitations are form a continuous spectrum. The pure Yang-Mills 
theories may belong to this class, according to the strong coupling expansion 
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[T3] . The result in these models is that ija/s ~ T 1+w where w is the power 
of the power law appearing in the density of states when expanded around the 
threshold. In generic case this ratio goes to zero at zero temperature. Therefore 
in these class of models there is no lower bound for the shear viscosity entropy 
density ratio. 
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